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Abstract 



The quantum of action h, multiplying in certain powers perturbative ver- 
tices in 4D gauge theory, is related to the action of just-not-resolved selfdual 
and thermal gauge field configurations, calorons and anticalorons, of charge 
modulus unity. Appealing to the derivation of the effective theory for the de- 
confining phase of SU(2) Yang-Mills thermodynamics, we conclude that these 
vertex inducers convey a rapidly decreasing interaction strength between fun- 
damental plane waves when the momentum transfer is increased away from 
the scale of maximal, effective resolution. This adds a deeper justification to 
the renormalization programme of perturbation theory which ignores the con- 
tribution to the partition function of nontrivially selfdual configurations. We 
also point out a connection between the QED fine-structure constant a and 
the electric-magnetically dual of the effective gauge coupling in the deconfin- 
ing phase, and we illustrate the workings of effective loops in the expansion 
of the pressure. 



1 Introduction 



Key experiments in all areas of subatomic and condensed matter physics, span- 
ning over more than a hundred years and designed to investigate aspects of the 
microscopic interaction between electromagnetic radiation, charged particles, and 
collective excitations, prove that the quantum of action h = 6.63 x 10~ 34 Js, intro- 
duced by Max Planck in understanding the spectral radiance of a black-body cavity 
through eventual appeal to the statistical thermodynamics of Ludwig Boltzmann pQ, 
is a universal constant of Nature. Since its inauguration in 1900 the somewhat mys- 
terious sneak-in of h, somehow being associated with 'elementary regions' or 'free 
rooms for action' [2], into the conception and development of quantum mechanics 
ID El IH El El E] has been broaching a quest for deeper interpretation. 

To understand the meaning of h as a mediator of indeterminism in fundamen- 
tal interactions among propagating gauge fields, between gauge and charged matter 
fields, or in stationary binding physics could be of relevance in addressing as of yet 
poorly understood phenomena, for example those of 2D strongly correlated elec- 
tron systems [8]. Quantum Mechanics and perturbatively accessed Quantum Field 
Theory are highly successful, accurate, and efficient. However, both suffer from 
problems of mere interpretation. While the measurement process in Quantum Me- 
chanics, that is, the interaction of a macroscopic observer with the subatomic system 
under investigation, mysteriously and in a nonlocal, instantaneous way collapses the 
evolved wave function of the latter to an eigenstate of the operator- valued, measured 
quantity the need for renormalization in Quantum Field Theory, in addition to not 
easing the classical Euler-Langrangian ambiguity inherent to the value of vacuum 
energy, is perceived by many as a necessary evil in putting perturbation theory to 
work. Fortunately, the predictive power of 4D gauge theories is not spoiled by the 
process of renormalization |HJ HU] , and radiatively introduced gauge anomalies can 
be canceled granting the consistency of the highly successful Standard Model of 
Particle Physics. Yet one may wonder whether the subtraction of infinities in the 
perturbative loop diagrams of gauge theories enjoys a deeper justification possibly 
supplied by full nonabelian gauge dynamics. 

In the present paper we consider the (Euclidean) action of a fundamental, self- 
duaQ gauge field configuration of topological charge modulus unity (trivial-holonomy 
Harrington-Shepard [TTj or nontrivial-holonomy Lee-Lu-Kraan- Van-Baal caloron/anti- 
caloron [12]) and scale parameter p in the deconfining phase of thermalized SU(2) 
Yang-Mills theory. The action of a caloron/anticaloron, that is just not resolved in 
the effective theory, is determined by the value of an effective coupling calculable 
in the according effective theory [T31 [TH EES] • Recall that calorons and anticalorons 
are not resolved in this effective theory: Their contribution to effective dynamics 
fundamentally separates into a free part causing the emergence of an effective, inert, 
and spatially homogeneous adjoint scalar field <fi [13J and the effects of interactions 
with the propagating plane- wave sector (vanishing topological charge). It is a fact 

1 Here the term 'selfdual' stands for both negative and positive topological charge. 
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that the energy-momentum tensor vanishes on any selfdual gauge-field configuration, 
see for example [TH], and hence also on calorons and anticalorons. Since the latter 
configurations thus do not propagate one may wonder [25] what their (potentially 
important) contribution to the thermodynamics of the high-temperature Yang-Mills 
gas actually is. 

For the generic nonabelian gauge group SU(2) this question was addressed in 
[T3" l [Tl j [To], see also [TB] for a summarizing presentation: While the emergence of 
is sharply dominated by (large) calorons/anticalorons of radius p smaller than but 
comparable^) to the scale of maximal resolution \(f>\ in the effective theory [13] inter- 
actions between calorons/anticalorons and plane waves manifest in a two-fold way 
effectively. First, there are small holonomy changes, induced through momentum 
transfers by plane-wave fluctuations that are largely off their free mass shell, which 
generate negative ground-state pressure by virtue of a pure-gauge configuration of 
the effective gauge field [HI [161 E] • Second, as we shall make explicit in this paper, 
just-not-resolved calorons/anticalorons are responsible for pointlike vertices, asso- 
ciated with their dissociation into screened but unresolved monopole-antimonopole 
pairs on one hand [IB] and a rapidly converging loop expansion in the effective the- 
ory [JJJ1 EDI Ell E2] ■ We will make clear that this loop expansion, which can naively 
be organized into an expansion in powers of h = turns out to be a more sub- 
tle expansion. Still, the action of a just-not-resolved caloron or anticaloron, which 
introduces a vertex for effective plane waves, turns out to be h. 

Calculating loops in effective variables yields a quantitative description of collec- 
tive effects as induced by fundamental field configurations [HI EHJ EH [25]. Moreover, 
since the fields of the effective theory and their interactions do only depend on the 
contributions of large calorons/anticalorons we would conclude that, on a funda- 
mental level, the mediation of interactions between unresolved off-shell plane-wave 
fluctuations by calorons/anticalorons (vertices) dies off rapidly with decreasing p. 
As a consequence, there are no interactions with potentially large momentum trans- 
fers deep within the unresolved physics. We tend to interpret this as a justification 
for the perturbative renormalization procedure in nonabelian gauge theories. Fi- 
nally, it is instructive to illustrate the technicalities of the effective loop expansion 
to contrast it with renormalized perturbation theory. 

The paper is organized as follows. Because our subsequent discussion relies on 
it we review in Sec. 12.11 the emergence of an inert, adjoint scalar field by a spa- 
tial coarse-graining over noninteracting calorons and anticalorons of charge modulus 
unity within the deconfining phase of SU(2) Yang-Mills thermodynamics. For the 
same reason, we re-derive in Sec. [272] the thermodynamically consistent solution for 
the temperature dependence of the effective gauge coupling e. In Sec. [273] we identify 
the Euclidean action of just- not-resolved calorons and anticalorons with the quan- 

2 The integral over p in the space-averaged two-point function of the fundamental field strength 
defining in singular gauge the kernel JC of a differential operator D, with <f> £ IC, exhibits a cubic 
dependence on the cutoff when the latter is comparable to |</>| -1 , see [13] and Sec. l2.1l below. 

3 All unresolved plane- wave fluctuations are at least by |0| off their free mass shell. 
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turn of action h and deduce, by virtue of the dimensionlessness of the fine-structure 
constant a, an interpretation of the fundamental unit of electric charge in terms of 
the dual to the effective SU(2) gauge coupling in the deconfining phase. Sec. 12.41 
discusses the implications of the action of a just-not-resolved caloron or anticaloron 
being associated with a local vertex. Namely, from the derivation of the effective 
theory we learn that such vertices are sharply suppressed with increasing momen- 
tum transfer (or decreasing caloron/anticaloron radius p) yielding, at any externally 
prescribed resolution, finite answers for all radiative processes. Focusing on the ex- 
ample of QED, Sec. 12. 51 discusses how just-not-resolved calorons and anticalorons of 
SU(2) cmb may induce local vertices between the plane waves belonging to electri- 
cally charged particles and propagating electromagnetic waves. To point out the 
contrast to renormalized perturbation theory, where the emergence of a scalar field 
associated with the maximal resolution and induced by the topologically nontrivial 
sector can not occur due to an a priori restriction to the plane-wave sector, Sec. [3] 
discusses the technicalities of the effective loop expansion of the pressure. Namely, 
we review the Feynman rules and constraints on loop variables and apply them to 
the (numerical) computation of three- and two-loop contributions to the pressure. 



2 Selfduality, propagating fields, and powers of h 
2.1 Inert adjoint scalar field and effective action 



In [13] the following definition for the kernel K of the differential operator D as- 
sociated with 0's equation of motion in the deconfining phase of SU(2) Yang-Mills 
thermodynamics was shown to be unique: 



K EE {0 a } EE J^tr / d 3 X ! 

C,A J J 



dpt a F^(r,0) {(r,0),(r,x)} F^(r,x) {(r, x), (r, 0)} , 

(1) 



where 



{(r,0),(r,x)}=Pexp 
{(t,x),(t,0)} = {(t,0),(t,x)? 



(r,x) 
(r,0) 



(2) 



and F^ v denotes the fundamental field-strength tensor. The Wilson lines in Eq. (j2J) 
are calculated along the straight spatial line connecting the points (r, 0) and (r, x), 
and V demands path-ordering symbol. (The euclidean time coordinate is r.) In 
(PQ) the sum is over a noninteracting caloron (C) and a noninteracting anticaloron 
(A) of trivial holonomy and topological charge modulus unity, and p denotes their 
(instanton) scale parameter. Both caloron and anticaloron are spatially centered at 
z c,a = 0, for a discussion of the distribution of action density in dependence of p/ ft 
(ft ee 1/T, T temperature) see [27]. In evaluating the integrals on the right-hand 
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side of Eq. ([T]) , one notices that only the magnetic-magnetic correlation gives rise to 
a nonvanishing result. Going over to spherical coordinates, both the unconstrained 
integration over r = |x| and p generate diverging multiplicative constants in front 
of a rapidly saturating sinusoidal time dependence separately arising from both the 
caloron as well as the anticaloron contribution. On the other hand, the azimuthal 
angular integration yields a muliplicative zero [T3~| IT6] . Together this generates an 
undetermined real prefactor which is one of the parameters spanning the kernel K, of 
a uniquely determined second-order linear operator D. More specifically, one arrives 
at the following expression for the kernel K. = {4> a } of differential operator D 

{0 a } = {Z c (5 al cosa c + 5 a2 sin a c ) A (27r(f + f c )) 

+E A (6 al cosa A + 6 a2 sma A )A(2n(r + f A ))} (a = 1,2, 3), (3) 

where He, E A , fc, and f A are undetermined real parameters, and ac, a A are angles 
that reflect the freedom of a global gauge choice. The amplitude function A is given 

as 

am ~\ 327r? fV~4 r r • /-«w- - -^1 vr 2 p 2 + cos(27rf) + 2 
A(2nr) = — - dpp 4 hmsm(2 5 r,r,p) — ) ' (4) 

3 Jo Lr^oo J (27T 2 /) 2 - COs(27Tr) + 1) 

where £ represents a cutoff for the p integration in units of (3. (All quantities of 
dimension length are made dimensionless by division with (3 and are denoted by 
a hat symbol). Function g saturates rapidly to a finite limit as p, f — > oo, the 
integral over p in Eq. (@| diverges cubically for £ — > oo, and function ^4(27rf) rapidly 
approaches constoo x £ 3 sin(27rf). Numerically, one obtains constoo = 272.018. To 



summarize, the kernel /C, which uniquely determines operator D as D = (9 2 + 

is sharply dominated by contributions of the p integration that lie in the vicinity of 
the upper limit £. 

Operator D is explicitly dependent on temperature and thus, by itself, cannot 
be taken to be associated with a temperature-independent effective actioro for the 
field <p. The latter, however, arises by allowing a potential V to absorb the ex- 
plicit temperature dependence introduced by D. Appealing to the BPS condition 
inherited from the selfduality of the fundamental gauge-field configurations entering 
into the definitition ([1]), a first-order equation for <fi is derived whose solution is 
parametrized by one constant of integration - the Yang-Mills scale A. Perturbative 
renormalizability [HI EI] , constraining the interactions between effective plane- wave 
fields a M = a^ta (t a with a = 1,2,3 the generators of SU(2) in fundamental repre- 
sentation, normalized to tr t a t b = \8 ab ), gauge invariance, and the fact that the field 
cannot, by definition, absorb or emit any energy-momentum yield the following 
unique result for the effective Euclidean action density £ cff [13l [16] 

(\ A 6 \ 
£ cff M = tr I- G^G^ + (D^) 2 + - , (5) 



4 This action is not explicitly temperature dependent because the weights of selfdual configura- 
tions in the fundamental partition function are solely associated with their topological charges. 
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where G^ u = d^a u — d v a^— ie[a^ : a u ] = G a ^ u t a denotes the field strength for gauge field 
a M , D^cj) = d^(p — ie[a^, 0], and e is the effective gauge coupling to be determined, see 
Sec. 12.21 Action density £ cff yields a highly accurate tree- level ground-state estimate 
and, as easily deduced in unitary gauge = 2|0|t 3 , tree-level mass m = 2e\<p\ = 



2 y ^ or P ro P a g a ti n g gauge modes a*' 2 rendering them thermal quasiparticles. 
2.2 Thermodynamical consistency 

For action density C cB to be thermodynamically consistent implied quantities like 
pressure and energy density must be related by Legendre transformations. On the 
level of free thermal fluctuations^] we obtain the following evolution equation 



24A 4 a D(2a) 
° (2tt)« 1 + ^fD(2a) ■ [b) 



Eq. (EJ) is equivalent to [26] 



24A 3 f . da i r , , 

A— + a aD(2a) , (7) 



(2tt) 6 V d\ 

where A = 2 4^-, a = ^ and thus a = 27reA~ 3//2 , and = f n °° c?x , x2 — , 9 1 , — . 

A ' 2T ' Jo y/x*+^ e Vx 2 +v 2 -l 

The right-hand side of Eq. (jSJ) is negative definite. As a consequence, its solution 
A (a) is strictly monotonic decreasing and so is its inverse a(A). Hence there is a 
regime in temperature A > Ai where a(A) < 1, and Eq. (JJJ) can be approximated as 

A 3 / \da \ , , 
1 = --— - A— + a o. (8) 



(2tt) 4 V 

The solution to Eq. (JHJ), subject to the initial condition a (A*) = Oj <C 1, is given as 

1/2 



^» = 4 ^ 3/2 ( 1 -^[ 1 -S) 



(9) 



Thus for Ai < A Aj function a(A) runs into the attractor a(A) = 4a/27t 2 A~ 3//2 . 
Using a = ^ = 27reA~ 3//2 , this attractor corresponds to the plateau e = y/8ir. 
Because the attractor increases with decreasing A the condition a <C 1 will be 
violated at small temperatures. The estimate 14.61 > Ai is obtained by setting the 
attractor equal to unity. Lowering A further, the solution to the full Eq. (J7J) runs 
into a thin pole at A c of the form [16] 



e(A) = -4.59 log(A - A c ) + 18.42 , (10) 



5 In the physical unitary-Coulomb gauge <f> — 2\<t>\t$ , 9,af = the off-shellness of quantum 
fluctuations in the field a M is constrained by the maximal resolution \<p\. Their contribution turns 
out to be negligible compared to that of thermal fluctuations O [16] . 
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where A c = 13.87. Now, setting the upper integration limit £ in Eq. (j3J) equal 
to l/(/3|0| -1 ) and invoking A c = 13.87, yields £ > (13.87) 3 / 2 /(2vr) ~ 8.22. But 
the approach of function A(2nf) to a sine is perfect for £ ~ 8.22 [T3J 02], and for 
£ ~ 8.22 the integral in Eq. (j3J) is strongly dominated by its upper limit. We conclude 
that calorons and anticalorons that are just not resolved dominate the emergence 
of the field <fi and its consequences for the effective gauge dynamics (adjoint Higgs 
mechanism and thermal ground state estimate). 

The divergence of e at A c initiates a transition to a phase where thermal quasi- 
particles decouple, the ground state is a condensate of magnetic monopoles, and a 
formerly massless gauge mode acquires mass by the dual Meissner effect [TH [16]. 

We conclude that almost everywhere in the deconfining phase the effective cou- 
pling assumes the constant value e = y/8ir. 



2.3 Counting of powers of h and action of just-not-resolved 
caloron 

Concerning the counting of powers in H, loop expansions in the effective theory 
are carried out as in conventional perturbation theory. Additional subtleties arise 
because of the existence of a maximal resolution \<p\ which enforces constraints on 
momentum transfers in four-vertices and on the off-shellness of the massless mode. 
These constraints have an obvious form in physical unitary- Coulomb gauge, see 
Sec.[3j To make the power counting in h most explicit, we work in units where 
kB = c = eo = fio = 1 but h is re-instated as an action. (So far we have worked in 
supernatural units h — c — ks — !)• The (dimensionless) exponential 



Iodrd 3 xC'J 



in the weight belonging to fluctuating fields in the partition function^] thus can, in 
unitary gauge, be re-cast as 

— J drd 3 xtr {^{d^a v - dyCb^ - iey/hla^ a u ]) 2 - e 2 h[a fl , 0] 2 j , (12) 

where d^ = a^/ y/h and (f) = (j)/\/h are assumed to not depend on H [29] , see also 
[30] . Notice that because of the terms oc h° in (1121) the unit of a M is length -1 . Thus 
the coupling e must have the unit of 1/y/h (and the unit of (j) is a ls° length -1 ). 
Together with the results reviewed in Sec. 12.21 we thus have 

(13) 



6 Recall that 4> is inert. As a consequence, the factor, whose exponent is the potential-part of 
the effective action, can be pulled out of the partition function and needs not be considered in a 
discussion of the effective loop expansion. 



6 



almost everywhere in the deconfining phase. Because fundamental field configura- 
tions are smoothly connected to the effective dynamics^ the action S^a-^ui- 1 °f a 
just-not-resolved caloron/anticaloron reads 

8tt 2 

Sc^p-m- 1 = -p- = ft- ( 14 ) 

Equation fJT^jl has implications. First, it suggests that the quantum of action h, 
whose introduction in flTTj) is motivated by the laws of Quantum Mechanics and 
which should enable a systematic accounting of quantum corrections (number of 
loops) in the effective theory, coincides with the Euclidean action of a just-not- 
resolved selfdual (classical) field configuration whose topological charge is fixed to 
be of modulus unity [13j IB]- That is, the reason why h really is a constant can 
be traced to the universal constancy of e (almost no dependence on T and no 
dependence on the specific Yang-Mills scale A of a given SU(2) gauge theory) and the 
fact that only charge-modulus-unity calorons/anticalorons contribute to the effective 
thermodynamical consistency of interacting, fundamental topological and plane- 
wave configurations [T3]. Moreover, rather than multiplying first and second powers 
of h 1 / 2 onto the powers e and e 2 of an a priori unknown coupling constant e in 
three- and four-vertices of a Yang-Mills theory we would infer that these vertices 
exist only because of the presence of just-not-resolved calorons and anticalorons. 
Second, for the fine-structure constant a of Quantum Electrodynamics (QED) to be 
dimensionless, 

'-^m- (18) 

the coupling g in Eq. ( Tl5|) must have the unit of yU. (N^ 1 denotes a numerical factor 
related to the mixing of the massless modes belonging to several SU(2) groups, see 
[16j|3T]-) This is guaranteed if g is taken to be the electric-magnetically dual coupling 
to e: 

4:71 r- 

g = — Kvfi. (16) 

e 

That is, a magnetic monopole liberated by the dissociation of a large-holonomy 
SU(2) caloron and other incarnations of this magnetic charge in the preconfining 
and confining phasd^] are interpreted as electric charges in the real world. Third, by 
virtue of Eq. (I13p the factors eyh and e 2 h in ( Tl2|) do not depend on h. Thus, by 
assumption the weight of the partition function associated with the effective theory 
is independent of h, and thus the effective loop expansion is not an expansion in 
powers of h. This is in contrast to renormalized perturbation theory where the value 
of the coupling is a priori unknown: The perturbative coupling constant needs to 
be determined by an empirical boundary condition. 



7 The spatial coarse-graining rapidly saturates: ^ > 8.221 x an d practically no 

alteration of 0's time dependence occurs when varying r, p ~ |16j . 

8 For example, the isolated charge situated in the center of the flux eddy associated with the 
selfintersection of a center- vortex loop in the confining phase Hi ,32^, see also discussion in Sec. 12.51 
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2.4 Unresolved, interacting plane waves: Local vertices by 
caloron/anticaloron mediation 



From the results of Sec. 12.31 and the facts reviewed in SecfU we may conclude that 
in the effective theory the pure-gauge configuration, associated with the a priori 
estimate of the thermal ground state, describes the collective effect of interacting 
(typically small-holonomy) calorons/anticalorons of action fi that are just not re- 
solved. Recall that the holonomy of calorons/anticalorons is elevated from small 
to large by effective plane-wave interactions which collectively account for dissoci- 
ation into screened but unresolved magnetic monopole-antimonopole pairs [T8"t 155] . 
Extrapolating the finding, that it is unresolved calorons/anticalorons which induce 
local vertices between plane waves, deeply into the unresolved domain, we would 
by the (ultraviolet as well as infrared) finiteness of all quantities calculable in the 
effective theory and the fact that the emergence of the field <ft is dominated by just- 
not-resolved calorons/anticalorons reason (cubic dependence of the prefactor of 0's 
time dependence on the cutoff for the p integration, see Eq. (jl]) and [T31 HE]) that 
the strength of interactions between fundamental plane waves, which are largely off 
their free mass shell, rapidly ceases with increasing momentum transfer. Recall that 
the effective pure-gauge configuration of the thermal ground-state estimate, which 
collectively describes these interactions, is sourced by the field [HI |T51 [IE]. But 
again, the content of fundamental charge-modulus unity selfdual gauge-field config- 
uration in (f) is strongly dominated by those calorons/anticalorons that are just not 
resolved. 

The point of view, see Fig.[H that the emergence of local vertices between plane 
waves in a Quantum Yang-Mills theory is a consequence of the existence of nonprop- 
agating, topologically stabilized and selfdual Euclidean field configurations of the 
same theory, compares in an interesting way to pure perturbation theory where such 
local vertices are taken for granted 0. Namely, while perturbation theory requires 
a renormalization programme to sweep ultraviolet divergences under the carpet the 
latter do not occur to begin with when vertices are understood to be induced by 
unresolved, nonpropagating, topological, and selfdual field configurations: Because 
potential plane-wave fluctuations, that are extremely off-shell compared to an ex- 
ternally imposed resolution^!, do not interact they are not generated in the first 
place. To the authors' minds this situation does by no means contradict the renor- 
malization programme of perturbation theory but rather supplies it with a deeper 
justification: The subtraction of infinities in perturbative loop diagrams, which ig- 
nore the presence of the nontrivially selfdual sector of gauge-field configurations, 
simulates the nonoccurrence of such divergences in the full Yang-Mills theory. 

9 For the effective and local interaction between fundamentally charged chiral fermions via their 
zero modes, which are localized on instantons, the 't Hooft vertex represents an example of plane 
waves interacting by the mediation of selfdual gauge-field configurations [34] [35] . 

10 In the case of thermodynamics this resolution is prescribed by temperature to be \<j>\. 
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(a) 



(b) 



Figure 1: (a) Plane waves (lines) of a Yang-Mills theory interacting via local vertices 
(solid circles), (b) generation of a vertex (local crumble within the ensemble of non- 
interacting plane waves) by a nonpropagating and selfdual gauge-field configuration 
of topological charge modulus unity in the same Yang-Mills theory. 

2.5 Gauge-field-matter interactions 

So far we have discussed the situation of a pure SU(2) Yang-Mills theory. The 
undisputed successes of perturbatively accessed QED in particular (and the Stan- 
dard Model of Particle Physics as a whole) demand to put the results of Sees. 12. 31 
and 12 .41 into the perspective of electrically charged particle interacting with photons. 
Let us undertake a perhaps speculative discussion of the situation in QED. 

To attribute caloron/anticaloron mediation to a QED vertex, the U(l) gauge 
group must be considered a dynamically broken SU(2). As explained in Sec. [Hand 
investigated in [131 HH HHJ US] this gauge-symmetry breaking is collectively mediated 
by calorons and anticalorons in terms of the field 0. In [36] the critical temperature 
T c of this SU(2) theory (and thus its Yang-Mills scale) was, by interpretation of the 
observed low-frequency excess in line temperature (see references in [36J), extracted 
to be the present baseline temperature 2.725 Kelvin of the Cosmic Microwave Back- 
ground (CMB). Thus the name SU(2) CMB . 

Radiative corrections in QED are organized into a power series in the fine- 
structure constant a. As already mentioned in Sec. 12. 31 a measures the strength 
of interaction of the massless mode of SU(2) CMB (the propagating photon) with a 
screened yet unresolved monopole [18J. Such a monopole becomes isolated and thus 
resolvable when trapped in the core of the selfintersection of a center-vortex loop 
[37] . Notice that in the presence of propagating photons a center- vortex loop of 
selfintersection number unity is the only absolutely stable degree of freedom in the 
confining phase of an SU(2) Yang-Mills theory. (The nonselfintersecting species 
shrinks to a round point [32] and vortex loops of intersection number larger than 
unity decay by charge annihilation or repulsion [371 [38] .) For the SU(2) CMB photon 
to see the charge of a monopole belonging to another SU(2) theory, which presently 
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is in its confining phase, the two theories must mix. So the crumble in the photonic 
plane wave, induced by the presence of the electric charge of, say, an electron, oc- 
curs by having the latter couple to its very massive photon (confining phase) which 
in turn mixes into the SU(2) CMB plane wave by caloron/anticaloron mediation of 
the latter theory. Recall that, according to the results of Sec. l2.3[ this mediation is 
accompanied by a unit h of 'free rooms for action'. 

Interestingly, since the QED action, involving the electron field ip and a^, is given 
as [2H] 

S QED = J d 4 x {^{ihd, -gal)-m)^, (17) 
where now ip , m = y , and are considered independent of h, one arrives at 

= (W(^ M - ±=.g oj) - roV (18) 

By virtue of Eq. ( |16p we would thus conclude from Eq. (ITS]) that the interaction 
■^gipj^a^ip does not exhibit an h dependence. So even in QED the loop expansion 
can no longer be considered an expansion in powers of h if the vertex is interpreted 
to be induced by a caloron/anticaloron of SU(2) CMB . 

3 Effective and irreducible, planar 3-loop correc- 
tions 

After 'scratching into unresolved physics' in Sec. [2] to associate interactions between 
plane waves with the presence of calorons/anticalorons we would here like to re- 
iterate [221 [39], HO] on technicalities in performing the effective loop expansion of the 
pressure in deconfining SU(2) Yang-Mills thermodynamics ]2T] . 

3.1 Feynman rules and constraints 

Following the extraction of the effective gauge coupling e(A) from the thermody- 
namical consistency of effective, noninteracting quasiparticles, one may investigate 
the magnitude of radiative effects implied by the so-obtained function e(A). This 
strategy turns out to be self consistent. To perform the calculation we go back to 
supernatural units and work within the real-time formulation [TT1 [T5] ITS] . Thus 
<V — > dfj, etc. 

The Feynman rules are formulated in the completely fixed (and physical) unitary- 
Coulomb gauge which avoids the consideration of Faddeev-Popov ghosts. 
• For the three- vertex rf^ 6c (p, k, q) and for the four- vertex r^j^ ,(p, q, r, s), we read 
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off from Eq. flTJJ that 

[3] abc 



T^(p,k,q) = e(27r) 4 8(p + q + k)e abc [g^(q-py + g^(k-qr + 



9^( P -kn (19) 



r%&Jp,Q,r,8) = -te 2 W5(p + q + s + r)ie fab e fcd (g» p g ua -g" a g up ) + 



efaceMg^g^ - g^g p n + e fad s fbc (g^g^ ~ g w g av )\ , 

(20) 

where e abc are the structure constants of SU(2). 

• The propagator of the tree-level heavy (TLH) modes propagates thermalized on- 
shell particles (no energy-momentum tranfer to <fi in the summation of the Dyson 
series for the interaction of a*' 2 with 0) of mass m = 2e\(f)\ only [16] 

>°l 

D„ = (^"^)- (22) 

The propagation of tree-level massless (TLM) modes splits into a vacuum and a 
thermal part: 



D™ b {p) = -2n6 ab D,J(p 2 -m 2 )n B ^ , (21) 



p 2 yjr ' " \ T 



- 2 + 2*6tf)n B [ ^ 



P 2 



(23) 



where 



poo = = = 0, P% = 6 ij - ^ . (24) 

Here ub(x) = denotes the Bose-Einstein distribution function. Notice that 
TLM modes have a color index a = 3 and TLH modes a = 1 or a = 2. Moreover, 
the term oc u^Uy in Eq. fl23l expresses the instantaneous "propagation" of the a% 
field (Coulomb term0), and = (1, 0, 0, 0) represents the four- velocity of the heat 
bath. 

• The maximal off-shellness of a TLM mode with four-momentum p is constrained 
as [HEE1CE6]: 

\p 2 \ < \(j)\ 2 (TLH modes) \p 2 \ < |0| 2 (TLM modes) . (25) 

For a four- vertex with ingoing four-momenta pi, pi and outgoing four- momenta P3,Pa 
we have [TBI EH EE] 



IG°i+P2) 2 | < \(j)\ 2 (s-channel), \(p 3 - pi) 2 \ < \(f)\ 2 (t-channel), 
\{Ps-P2) 2 \ < \<P\ 2 (u-channel). (26) 



11 However, one-loop radiatively induced, the emergence of longitudinally propagating £?-field 
(charge-density) waves is observed [24] . 
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As we shall review [21] [22] in Sec. l3.2[ conditions ( 125]) . ( 126|) . and the requirement of 
on-shellness of TLH modes imply a rapid growth of the number of constraints over 
the number of a priori noncompact (radial) integration variables with an increasing 
number of loops in an irreducible diagram. 

3.2 Expectations for convergence properties of loop expan- 
sion from two- and three-loop corrections to the pres- 
sure 

As discussed in the previous section, the momentum transfer in all Mandelstam 
variables s,t, and u in a four-vertex with ingoing four-momenta pi,p2 and outgoing 
four-momenta p^,P4 should not exceed the scale \<f>\. For three- loop diagrams, and 
in particular the ones of Fig.[3j where the cutting of any line yields a lPi irreducible 
contibution to a polarization tensor, the number K of independent, radial loop- 
momentum variables (po, \p\)i for % = 1, 2, 3 is K = 6. To judge the number K 
of independent constraints imposed on them we, in addition to (|25|) and (|26|) . to 
consider on-shellness constraints for the TLH modes: 

p\ = m 2 , p\ = m 2 , p\ = m 2 , p\ = (p 1 + p 2 - p 3 ) 2 = m 2 . (27) 

Conditions (1261) and (I2"7j) together yield for irreducible three-loop diagrams A and 
B a total ofi^ = 3 + 4 = 7 constraints. 

According to (125]) . the TLM modes in the irreducible three- loop diagram C are 
subject to the maximal off-shellness conditions 

b?|<|0| 2 , \PI\<\^\ 2 , (28) 
and the on-shellness relations for the TLH modes in diagram C are 

p 2 3 = m 2 , p\ = (px +p 2 - P3) 2 = m 2 . (29) 

Thus, again, we have K = 3 + 4 = 7 constraints for diagram C. Hence for all 
irreducible three-loop diagrams we have K = 7. Therefore, the number of constraints 
exceeds the number of radial variables: 

K<K. (30) 

This relation implies that irreducible three-loop integrations have either compact or 
empty supports. 

For the two-loop diagrams of Fig. [2], we have K > K. Namely, the number of 
K of independent radial loop momenta (po, \p\)i for i = 1,2 is K = A whereas the 
constraints plus the on-shellness conditions are 

\(pi+P2) 2 \ < |0| 2 , Pl = m 2 , p 2 2 =m 2 . (31) 
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P3 = P! + P 2 

A B 

Figure 2: Two out of three nonvanishing two-loop diagrams with TLH (solid) and 
TLM (dashed) modes in the loops. 




P4 = Pj+ P2T P3 P4 = P|+ P2" ?3 ?4 = Pj+ P2" ?3 

P5=P r P 3 

A B C 



Figure 3: The three irreducible, planar three- loop diagrams with TLH (solid) and 
TLM (dashed) modes in the loops. 
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This implies K — 1 + 2 = 3. Thus, unlike the case of irreducible three- loop integra- 
tions, we expect the support of two- loop integrations to be noncompact. In Ref. [2TJ 
a counting of K versus K was performed for an arbitrary loop number in the case 
of planar and nonplanar irreducible loop diagrams which gives rise to the conjecture 
that, modulo resummations of irreducible contributions to the polarization tensors 
in the lines of a given bubble diagram, the expansion of the pressure terminates at 
a finite loop order. The fact that K < K [22] is already in agreement with this 
conjecture. 

In the remaining subsections of Sec. [3] we discuss how numerical results are ob- 
tained and that, compared to the free quasiparticle pressure, two-loop and irre- 
ducible three-loop diagrams exhibit a hierarchical suppression of their contributions 
to the thermodynamical pressure [22]- In fact, diagram C of Fig. [3] is shown to 
vanish. 

3.3 Irreducible 3-loop integrations 

For the irreducible three- loop bubble diagrams in Fig. [31 we have in case A two 
four-vertices and four propagators and the correction to the pressure is 

If d i pid A p 2 d A p 3 
LXt ^ A ~ 48 J (27r) 4 (27r) 4 )(27r) 4 ^ abcd \^bcd u PP^KPijV^^AV'i) y 6A ) 

where the momenta are subject to the constraints (1261) and (T2~Tj) and overall momen- 
tum conservation has been used (p 4 is a function of pi, p 2 , p 3 .). 

For diagram B we have one four-vertex, two three-vertices and five propagators. 
Therefore, the correction to the pressure is 

A p 1 f d pid p 2 d p 3 rQ/ 3 7 A 

^ rB ~ 2 J (27r) 4 (27r) 4 )(27r) 4 W hi i k [ 3 l abc Wabc^^h^^^MPz) {^) 

X D ~1P J c (P3 ) D X p, kc (P4 ) D^ D <h l (p 5 ) , 

where the momenta are subject to same constrains as in diagram A. It should be 
noticed that for the massless mode, p 5 = p\ — p 3 , we do not need to consider an 
additional constraint |]?5| 2 < |0| 2 because this constraint is automatically satisfied 
by the t-channel condition in ( J26l) . 

For diagram C, the Feynman rules are similar to diagram A and the correction 
to the pressure is 

AP 1 f d± Pl d%d 4 p 3 , m 

48 J (27r) 4 (27r) 4 )(27r) 4 W a6cd pAc5 s,<id 

x D^aaips) D uPtb i (p 4 ) , (34) 

where the momenta are subject to the constraints ( 1261) . ( I28i) . and ( 1291) . 
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Now, insert the Feynman rules (jT9l) - (1231) in the above integrals (1321) - (1341) . 
consider the dimensionless variables Xj = W (i = 1,2,3) for rescaling the radial 
loop momenta, and consider the following steps: 

• Carry out Lorenz and color contractions. 

• Go over to 3D spherical coordinates. 

• Trivially integrate the temporal components of four-momenta by appealing to the 
delta functions of the propagators. 

• Apply the triangle inequality to obtain upper bounds for the modulus of each 
correction. 

For the corrections (1321) and fl33l) we obtain the following expressions 



e 4 A 4 A~ 2 f f f f f r z 23,u 

APa(b) I < (24) 3 x 2 7 x (277)6 J dxi J dx2 J dx ' 3 J dzi2 J dzi3 J dz2 



''223,1 

Lm,n=l / J J J J Jz 23 ,i 



Vi 1 ' z h) (! - z h) ~ (^23 - z 12 z 13 ) 2 

111 

ty* fy HP 

X 1 X 2 X 3 

yjx\ + te 2 \/x\ + 4e 2 A/x| + 4e 2 
x5(4e 2 + {{-l) l+m ^/x 2 + Ae 2 ^/x 2 + Ae 2 - Xl x 2 z 12 - 

((-l) l+n ^x 2 + 4e 2 ^x 2 + 4e 2 - Xl x 3 z 13 ) - 

((-l) m+n V^ + 4 eV^i + 4e 2 - x 2 x 3 z 23 ) 

x |Pa(b)(^> ^ m -, n)\n B (27i\~ 3 / 2 ^ xj + 4e 2 x 

n B (27rA- 3/2 - v /x 2 + 4e% B (2vrA~ 3/2 -\/x 2 + 4e 2 x 

n B (27rA- 3 / 2 |(-l)yx 2 + 4e 2 + (-l) m ^/x 2 + 4e 2 + 

(-1)^x1 + 4e 2 |, (35) 

where Z\ 2 = cos Z(x 1; x 2 ), 243 = cos Z(x 1; x 3 ), 223 = cosZ(x 2 ,x 3 ). The extremely 
lengthy polynomials Pa and Pg arise from Lorenz and color contractions and are 
regular at x\ = x 2 = x 3 = 0. We also define 



z 23:U = cos |arccos2;i2 — arccoszi 3 | , 2:23^ = cos |arccos2;i2 + arccos^isl . (36) 



Recalling ( 127)) and the fact that only the minus sign is relevant (e > 
which makes the expression within the absolute-value signs strictly negative, we 
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have 



z 12 < -^-(4e 2 - ^x\ + Ae^x 2 2 + 4e 2 + i) = x 2 ) (37) 
z 13 > ^-(-4e 2 + ^xf + 4 e yx| + 4e 2 - i) = x 3 ) (38) 
z 23 > — (-4e 2 + 7x 2 + 4e 2 ,/x 2 + 4e 2 - ~) = p(x 2 , x 3 ) • (39) 

X 2 £ 3 V V i 

Repeating the same steps as before, for the vacuum-vacuum contribution f l34|) of 
diagram C we obtain 

2 



APc| < 3 x 2 5 x f27r1 8 ^ / rfX2 / dXs dZl2 l dZl3 dZ<23 

l,m "' "' z i'i,l 



111 

rf 44 nf" rrt* 4 

1 2 3 



x| 



|P c (x, z, yi , Z, m) |n B (2vrA- 3 / 2 Jx 2 3 + 4e 2 ) 



x ^(2 7 rA- 3 / 2 |(-iy v ^xfT4^+ (-ir/ 2 (x,z)|) ^ 
/2(x,z,yi)^/x§ + 4e 2 



where 



/ 2 (x,z) = y x 2 + x| + x| + 2xix 2 zi 2 - 2xix 3 zi 3 - 2a; 2 x 3 z 23 , 2/1 = ^ , (41) 

with 2!i2, -2i 3 , -223, ^23,«, and 2 23i ; defined in fl3T|) .f l3"8"]) .f l39~]) and fl36|) . respectively. The 
polynomial Pc emerges from Lorenz and color contractions and is regular at x\ = 
X2 = X3 = 0. According to relation ( )26l) - ( 128]) the rescaled constraints read: 





\yl + yl - 


■ 3C\ 


x\ + 2/12/2 - 2x 1 x 2 z 12 \ 


< 


1- 


(42) 


\y\ 


-x 2 + 4e 2 -(-i; 




Jx\ + 4e 2 + 2x 2 x 3 z 23 


< 


1- 


(43) 




-x 2 + 4e 2 - (-1) 




Jxl + 4e 2 + 2xix 3 ^i 3 


< 


1- 


(44) 






\vl- 


x 2 | < 1, \yl-x 2 2 \ 


< 


1. 


(45) 



Moreover, y 2 is replaced by 

2/2 = -2/1 + 2(-l)y* 2 + 4e 2 + (-ir/ 2 (x,z), (46) 
through delta function integration, and (l,m = 1,3). 
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3.4 Monte-Carlo integration of irreducible 3-loop diagrams 

3.4.1 Monte-Carlo integration in general 

In the following we briefly discuss how the integral of a function / over a region 
G C M. n can be calculated using statistical methods. Du to the complexity of 
the region (algebraic varieties) G that support^ a three-loop integration, statistical 
methods are required in a numerical computations of the values of bubble diagrams 
A, B, and C in Fig. [31 Its characteristic function is denoted by xg- The integral 
then reads 

/ f(x)dx = / xa(x)f(x)dx. (47) 

JG JR n 

If G is compact then it can be included in a box B = [xi,Xi] x [x2, X 2 ] X ... x [x n , X n ]. 
Therefore the integral in Eq. (T4"Tj) can be written as 

/ XG(x)f(x)dx = / xa{x)f{x)dx. (48) 

Jl" J B 

If B has volume V, ^Xb can be considered as the probability density function of a 
random variable X, which is equally distributed thoughout the box B. The integral 
then is exactly the expectation value E(VxG{X)f(X)) 

! X G{x)f{x)dx = V f XG(x)f(x)^ X Bdx = VE( XG (X)f(X)). (49) 

J B J B 

The Monte-Carlo method to calculate this integral consists of a statistical estimate 
of the expectation value. It is well known that the mean value X of a sample 
is an unbiased estimator for the expected value. If one draws a random sample 
(xi,X2---,x n ) of points from the box, the estimate becomes: 

f 1 n 

/ f(x)dx^V-y2xG(x i )f(x i ). (50) 
Jg njr( 

The Monte-Carlo-Method is particularly useful to determine integrals over high- 
dimensional integration regions, where deterministic methods would be too time 
consuming. Unfortunately, the precision of this statistical estimate increases only 
with the root of the sample size. For example, this implies that the Monte-Carlo 
method is a relatively fast way to achieve a result with a relative precision of about 
1 percent, but to achieve one more decimal place in the result the sample size must 
be increased by a factor of 100. 

3.4.2 Monte-Carlo integration for diagrams A and B 

The calculation of AP4 and APg includes an integration over a six-dimensional re- 
gion G with variables of integration El Xi, x 2 , x 3 , z i2 , Zi 3 , z 23 . As = cos Z(Xj, Xj), 

12 The region G is determined by a set of inequalities. 

13 Notice that x\ is integrated analytically through the delta function in APa and APb- 
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|AP a |/Pi_i 00 p 



6xl(T 14 
5xl0" 14 
4 x 10" 14 
3 x 10~ 14 
2 x 1(T 14 



1x10 



-14 




20 40 60 80 100 120 140 



Figure 4: Upper estimate for the modulus -^\APa\/ '-P]_i 00 p as a function of A for 
diagram A in Fig.EJ 



their range is restricted to the interval [—1,1]. The exact shape of G is then de- 
termined by the constraints (1371) - (1391) and an additional angular condition (|36|) 
on Z23. This angular condition allows for a further reduction of the range of Z23 to 
the interval [223,/, 223,«]- Also, the constraints ( |371) - (1391) and condition (1371) restrict 
the to the interval [0,3] [21] . Therefore a box I? which contains G is identified. 
The constraints define the characteristic function \g whereas z 2 3,« and z 2 3 )U are used 
directly as integration limits. According to this, the ratios of the moduli AP^ and 
APb to the one-loop (ground-state subtracted) approximation are computed as a 
function of the dimensionless temperature A c = 13.87 < A = < 140, and the 
results are represented in Fig.HHTJ For the effective gauge coupling the solution of 
the one- loop evolution equation ([7]) is used. Figs.0] - [7J show that at the critical 
temperature A c = 13.87 there is no loop correction (decoupling of TLH modes), and 
AP4 and APb approach zero for large temperatures. For the three-loop diagrams 
A and B, the maxima of the ratio of the estimates of their moduli to the one-loop 
correction are peaked between A = 17.5 and A = 20. They are 6 • 10~ 14 and 2 • 10 -7 , 
respectively. This shows the dominance of diagram B. (The next section suggests 
that diagram C is vanishing). Comparing Fig. |5] with Fig. |6] and Fig.0 shows the 
dominance of Coulomb fluctuations (10~ 7 ) over quantum fluctuation (10 -12 ) for in 
the TLM propagation contributing to diagram B. 



18 



|AP b l/Pi-i 00 p 




20 40 60 80 100 120 140 A 

Figure 5: Total upper estimate for the modulus ^\AP B \ / 'Pi_\ 00 p as a function of A 
for diagram B in Fig.[3j 

3.4.3 Monte-Carlo integration of diagram C 

The region of support for the radial loop integration the irreducible three-loop dia- 
gram C cannot be determined in the same way as for diagrams A and B. This has 
mainly to do with the occurrence of additional off-shell variables. In order to esti- 
mate the region of loop integration G for diagram C, we consider a box with volume 
10 x 10 x 10 x 20 x 2 x 2 x 2 for the loop and angular variables xi, x 2 , x 3 , yi, Zi 2 , Zi 3 , z 23 , 
respectively. Through the condition \y\ — x\\ < 1, and the definition of z 2 3 J ( u ,i) in 
Eq. (13 6 p a subregion of this volume is determined, which roughly represents 2 % of 
the volume of the box. In this subset 150,000,000 points are randomly picked. Sub- 
sequently, it is checked whether the conditions ( 142|) - ( 146]) . With 600,000,000 tests we 
have not found a single point to satisfy all these conditions. This yields an estimate 
of a fraction of about 2 x 10~ 5 of the box volume for the region of integration. The 
typical distance between the sampled point is 0.2 ~ \/2 x 10~ 5 . Thus it is highly 
probable that the region of integration is empty. A similar analysis applies to the 
vacuum-thermal combination of TLM propagators, and the constraints imply that 
the integration region is empty provided the region of integration for the vacuum- 
vacuum combination is empty. Again, we have gathered compelling evidence for the 
latter to be true. 

3.5 Comparison to 2-loops results, hierarchy 

In order to check the reliability of our numerical method in case of three-loop in- 
tegrations, we compute the pressure of the two- loop diagram B in Fig. [2] with the 
Monte-Carlo method and compare the result with results obtained by deterministic 
means [391 SO]. 
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I AP b |/Pi_ioop 

1 . 75 x 10" 7 
1 . 5 x 1(T 7 



1.25x10 
1x10 

7.5x10 
5x10 

2.5x10 




Figure 6: Upper estimate for the modulus of I APb |/-P]_i 00 p due to the Coulomb 
part of the TLM propagator as a function of A for diagram B. 



I AP b I /Pi-ioop 




20 40 60 



Figure 7: Upper estimate for the modulus of ^\APb\/ Pi_\ OQ p due to TLM quantum 
fluctuations function of A for diagram B. 
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For the two- loop diagram B in Fig.[2j we have two-propagators and one four- 
vertex at which the momenta are subject to the s-channel constraint constraint in 
(j2"Sj) . Thus the two-loop correction to the free quasiparticle pressure arising from 
diagram B takes the form 

1 f dtptfk^pg 



APb = -J -^ T SL D ^p)D P sAk) , (51) 

where 

|(p + fc) 2 |<|0| 2 . (52) 
As before, we can rescale the momenta as x = j^j and y = l^j and obtain 

e 2 A 4 A~ 2 r J r J r J xy 

AP B = — 0/r> _, 4 dx dy dz xy 



x 



2(2vr) 4 7 7 tf 7 Vx 2 + 4eV^ 2 + 4e 2 
P B (x, y, z xy )n B (2n\- 3 /Wx 2 + 4e 2 K(2vrA- 3/2 vV + 4e 2 ) , (53) 



where 



Notice that 



z* v < — (4e 2 - Vx 2 + 4e 2 vV + 4 e 2 + r) = g(x, y) . (54) 
xy 2 

lim g(x,y) = — 1. (55) 

In (|52|) only the minus sign is relevant since e > and therefore the expression 
within the absolute- value signs is strictly negative. This leads to ( 154"|) . Apart from 
a small compact region, where g(x,y) > 1 and which includes the origin x = y = 
in the (x > 0,y > 0)-quadrant, the admissible region of x, y-integration (—1 < 
g(x,y) — +1) is an infinite strip bounded by two functions: 



x + 8e 2 + y/1 + 16e 2 Vx 2 + 4e 2 
~8? 



yU{x) ^ -rv^ ^ V x ^ ) (56) 



X 



+ 8e 2 - VI + 16e 2 Vx 2 + 4e 2 



y'(x) = v ^ P v - ' - • (57) 

We now use the Monte-Carlo Method to calculate diagram B. In order to bound the 
region we restrict the arguments of the Bose-factors to the following interval: 



< 2tt\- 3/2 Vx 2 + 4e 2 < 10 (58) 

and 

< 27rA~ 3/2 vV + 4e 2 < 10 . (59) 

These restrictions ensure that no Boltzmann tails are included in the region of radial 
loop integration. The integrand of (I53p is plotted in Fig. [5] below. It is clear from 
Fig. M that the support for the x and y integrations is an infinite strip in the x-y 
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200 




Figure 8: The integrand in Eq. (1531) plotted as a function of x and y for fixed 
z xy = cos Z(x, y). The horizontal plane represents the x-y plane where the integrand 
(mountain-shaped object) stands on an infinite strip. 



plane. From x = y = 150 on suppressions are severe justifying the limits taken in 
(158]) and fl59|) . For instance, by taking x = y = 150, e = \/8it and A = 30, we obtain 
2%\~ 3 ^ 2 \/x 2 + 4e 2 ~ 5.77. The ratio of two-loop to one-loop, as computed by the 
Monte-Carlo method, is represented in Fig.[9j There are no significant deviations 
from the deterministic results computed in [SHI SO]- Both estimates are of order 10 -6 
at temperatures A = 20 ± 5. Let us now compare the integrand of (j5"3"|) . illustrated 
in Fig.[8]for the two-loop diagram B, with the integrand of (135]) . illustrated in Fig.fTUl 
for the three- loop diagram B. 

To obtain an idea about the shape of the region of radial loop integration for 
the three-loop diagram B we adopt the following procedure. We use a grid of cosine 
values for z 12 , zi 3 , and z 23 of 0.15 separation. This amounts to approximatively 2300 
points. For each point we compute the modulus of the integrand as a function of x 2 
and X3] the xi-integration can be performed trivially by appealing to the 5-function 
in the integrand. Subsequently, a sum over all grid points is performed, and Fig.fTOl 
depicts the result of this summation and the x\ integration. The x 2 -x 3 region of 
integration is the area where this result does not vanish. Fig. [TD] supports the claim 
in [2T] that the region of radial loop integration in the three- loop case is compact. 
The fact that the upper limits of integration for x 2 and x 3 have been set to 3 shows 
that this region is fully covered in the Monte-Carlo process. 
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Figure 9: The ratio of AP B (2-loop) and Pi_ ioop plotted for 13.87 < A < 140. 



Recall that Fig. [8] shows the integrand for the two-loop diagram for a particular 
choice of the cosine value z xy close to —1. The region of radial loop integration is 
obviously not compact. However, the integrand is Boltzmann suppressed for large 
x and y. 

The remarkable contrast between the two-loop and irreducible, planar three-loop 
diagrams can be seen from Fig. [8] and Fig.[10j in agreement with our counting of 
constraints in Sec. 3.2, the region of radial loop integration is compact in the latter 
while it is non-compact in the former case. The numerical results [4TH [22] , reviewed 
in Fig. H]- [7] and Fig. [9], show that one-loop, two-loop, and irreducible three-loop cor- 
rections are hierarchically ordered with a pronounced suppression at increasing loop 
order. In particular, comparing the modulus of the dominant irreducible three-loop 
correction with the smallest two-loop correction, we infer that they are comparable 
at their maxima in A. On the other hand, comparing it with the dominant two-loop 
correction, shows a minimal suppression by a factor of ~ 1CT 3 . The dominant irre- 
ducible three-loop contribution is significantly dominated by Coulomb fluctuations, 
and the other nonvanishing irreducible three-loop contribution is minimally sup- 
pressed by a factor 1CT 11 compared to the dominant two-loop contribution. Also, 
recall that one three-loop diagram appears to vanish exactly. These results pro- 
vide numerical evidence for the conjecture [21] on the termination of expansion of 
thermodynamical quantities or polarization tensors into irreducible loops (modulo 
resummation of reducible diagrams) to be true M . 



14 As discussed in [21j . constraints on the loop- integration variables in nonplanar, irreducible 
bubble diagrams arc more severe than in the here-discussed planar case. 
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Figure 10: The result of xi-integration and summation over a grid of cosine values 
z\2 = cosZ(xi,x 2 ), ^13 = cosZ(xi,x 3 ), and Z23 = cosZ(x 2 ,x 3 ) of the integrand in 
Eq. (ESD plotted function of x 2 and X3, see text. 



4 Summary 

In this work we have addressed the question how the universality of the quantum of 
action h may find an explanation in terms of the Euclidean action of nonpropagating, 
selfdual gauge-field configurations of an SU(2) Yang-Mills theory and, after inter- 
preting these configurations as the inducers of local vertices between plane waves, 
how the radiative correction of the full theory become finite. Recall that in per- 
turbation theory, where the topologically nontrivial sector of field configurations is 
ignored, ultraviolet divergences are removed by a conditioned subtraction. We have 
pointed out that this pragmatism of renormalization theory, which was proved a long 
time ago to be physically consistent [9j [10], is justified by solely a small window in 
scale parameter, associated with the just-not-resolved, relevant topological sector 
of the theory (calorons/anticalorons of topological charge modulus unity), inducing 
local vertices between plane waves. We also have discussed how the coupling of 
photons to electric charges is related to the dual of the effective gauge coupling in 
the deconfining phase. Finally, to illustrate the workings of effective loops we have 
made explicit certain technicalities of the effective pressure expansion in deconfin- 
ing SU(2) Yang-Mills thermodynamics. Numerically, the hierarchical suppression of 
increasing loop orders supports the conjecture of a termination at a finite order of 
the expansion into irreducible diagrams [2T] . 
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